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PREFACE 

This project is under the supervision of Professor Ivan I. Mueller, 
Department of Geodetic Science, OSU, and it is under the technical 
direction of Mr. Jerome D. Rosenberg, Deputy Director, Communications 
Programs, OSSA, NASA Headquarters, Washington, D. C. The contract 
is administered by the Office of University Affairs, NASA, Washington, 

D. C. 20546. 
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1. STATEMENT OF WORK 

The statement of work for this project includes data analysis and 
supporting research in connection with the following broad objectives: 

(1) Provide a precise and accurate geometric description of 
the earth's surface. 

(2) Provide a precise and accurate mathematical description 
of the earth's gravitational field. 

(3) Determine time variations of the geometry of the ocean 
surface, the solid earth, the gravity field, and other geo- 
physical parameters. 
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2. ACCOMPLISHMENTS DURING THE REPORT PERIOD 

2. 1 Adjustment of BC-4 Worldwide Geometric 
Satellite Triangulation Net 

During this reporting period, NOAA began the shipments to the 
Data Center of optical observations from the BC-4 worldwide network. 
This data is in two different forms and is referred to as Type I and 
Type II data. The Type I data gives satellite image plate coordinates 
for every satellite image on a plate. Along with this are the camera 
calibration parameters and a covariance matrix. The Type II data is 
the result of a polynomial fit to the plate images, which gives seven 
ficticious images on each plate. These ficticious images are selected 
at the same time for all stations observing the satellite so that for a 
particular event, the seven satellite images are simultaneous for all 
camera stations observing. The term "event", as used by NOAA, 
refers to a series of images on camera plates; there can be as many 
as four stations observing, and as many as 400 Type I images or 7 
Type II images on each plate. All of this data comprises one event. 

Due to the fact that all Type II data is simultaneous, NOAA used 
Greenwich Hour Angle and Declination as observations for each image. 

In this way, it was not necessary to record the time of observation. 
Because of the correlation among the ficticious images on a single 
plate, a 14 x 14 variance-covariance matrix is given for each set of 
plate data. 

The correlation between different satellite images cannot be taken 
into consideration in our existing optical adjustment program. In order 
to test the Type II data, we neglected the correlation between the images 
on a small sample of the data, generated a ficticious time and converted 
the GHA's to Right Ascensions, and then input this data into our existing 
adjustment program. The results showed that the observations 
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were not Greenwich Hour Angle and Declination. A telephone call to 
NOAA in Rockville, Maryland verified that the observations were not 
GHA and Declination, but actually Azimuth and Elevation with respect 
to the origin station, Aberdeen. NOAA then sent a letter describing 
how we could convert this data to GHA and 6. This conversion was 
made in our program and used with the old data tapes until NOAA 
actually sent the correct tapes. 

When the correct Type II data tapes finally became available (in 
December) we had already written and had operational a computer 
program that would read the data (without correlation) from the 
original tapes and punch this data (Greenwich Hour Angle and Declina- 
tion) on cards for input into our existing adjustment program. Changes 
were also made in the existing adjustment program to accept the 
Greenwich Hour Angle directly without converting to Right Ascension 
as required originally. If we neglect the correlation between images 
we could process all of the data immediately. 

Because of the correlation between the images on each BC-4 
plate, it became necessary to write a new computer program to process 
these observations. In the existing program the normal equations were 
formed one event at a time, but in that case an event was all stations 
observing one satellite point. It is now necessary to form the normal 
equations using seven satellite points at once, and this required major 
changes because now the seven images from each plate and their 
associated 14 X 14 variance-covariance matrix must be processed 
together. This necessitated changing the linearized form of the math 
model. All of the logic has been developed and most of the program- 
ming is completed at this time. 


Attachments (2): 

Correspondence dated November 17, 1971 
and November 29, 1971 
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ATTACHMENT #1 



ATTACHMENT TO SECTION 2. It (6 Pages) 

U.S. DEPARTMENT OF COMMERCE 

National Oceanic and Atmospheric Administration 

Rockville, Md. 20852 

NATIONAL OCEAN SURVEY 



November 1 7, 1971 

Professor Ivan I. Mueller 
Dept, of Geodetic Science 
The Ohio State University 
164 West 19th Avenue 
Columbus, Ohio 43210 

Dear Professor Mueller: 

Dr. Schmid has probably told you, as you had already dis- 
covered, that the type II data we sent to NASA is not what 
it is supposed to be. It is supposed to be Greenwich hour 
angle and declination referred to the mean pole of 1900-1905, 
the Conventional International Origin. Just a word of ex- 
planation of how this came about, then I will give you the 
data and formulas needed to transform the data you now have 
into GHA and declination and also suggest bn alternative. 

When the former Coast and Geodetic Survey started experiment- 
ing with satellite triangulation in the early 1960's the 
first observations were made from a small triangle at Aberdeen 
Proving Grounds in Aberdeen, Maryland. The cartesian coordinate 
system established for triangulation had as its origin a 
point on the ellipsoid directly beneath one of the camera 
stations. The system was left-handed with the Z axis directed 
positive toward the zenith, the X axis positive toward north, 
and the Y axis positive toward east. The geodetic ellipsoidal 
coordinates of the stations were referred to the 1927 North 
American Datum, Clarke 1866 ellipsoid. 

As the experiments grew from small to large triangles, and 
even when projects continental and worldwide in scope were 
undertaken, the coordinate systems remained the same. 

In recent years the entire data reduction system was modified 
extensively and it was decided to reprocess all of the data. 

It was also decided that a new reference ellipsoid, viz. 

Navy 8-D, would be adopted and that an ellipsoid-centered 
cartesian coordinate system would henceforth be used. 
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At the time I wrote the programs to transform our data into 
type I and type II data for the NASA library, they were 
designed to be used exclusively with the reprocessed data, 
i.e. an ellipsoid-centered system. Having no reprocessed 
data available at the time, the programs were tested with 
fictitious data and functioned properly. 

As it turned out, for whatever reasons, the reprocessed data 
were referred to the Navy 8-D ellipsoid, as planned, but the 
cartesian coordinates continued to have as their origin the 
point at Aberdeen and an orientation similar to the earlier 
one (X north, etc.). Consequently, when these data were used 
to generate type II data for NASA, the quantities which we 
called GHA and declination were actually, or could be thought 
of as, an azimuth from north and -air elevation angle referred 
to the cartesian origin, the point^at Aberdeen. 

In order to salvage the existing erroneous type II data, 
one would have to refer the satellite directions to an 
ellipsoid-centered system by means of rotations and recompute 
the GHA-declination components. Details of this procedure, 
together with the Navy 8-D coordinates of the Aberdeen origin 
point are given in the enclosure for your inspection. 

Since your phone call of last Friday, I have taken steps to 
correct the error. The program for computing type II data 
has been modified to compute the proper quantities. All type 
II data previously sent to NASA will be recalled and will 
be replaced by new data. 

The alternative I mentioned earlier is for you to wait for 
correct data, rather than salvage the data you now have. 

Since all of the previous input to the type II data program 
has been recorded on magnetic tape, the correction process 
should be completed with' r 'a few weeks, at most. 

One other item which may also concern you has to do with 
type I data. "While there is nothing wrong, that we know of, 
with the data themselves, the portion of the associated tech- 
nical report describing the computation of GHA and declination 
from these data is, for the given reasons, also incorrect. 
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A revised report will be sent to NASA with a request that 
old copies be destroyed and users notified accordingly. 

Sincerely, 

I \ . 

Robert H. Hanson 
Geodetic Research and 
Development Laboratory 

Enclosure 
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ATTACHMENT TO SECTION 2.1: (15 Pages) 

U.S. DEPARTMENT OF COMMERCE 

National Oceanic and Atmospheric Administration 

Rockville, Md. 20852 

NATIONAL OCEAN SURVEY 


ATTACHMENT #2 


November 29, 1971 

Hr. Jerome D. Rosenberg 
Deputy Director, Communications Programs 
Office of Space Science and Applications 
National Aeronautics and Space Administration 
Washington, D. C. 205^6 

Dear Mr. Rosenberg: 

Attached is a revised version of ’the technical , report for 
'Contract W-13,321, dated January 19, 1971* Please note 
that only those pages containing changes have been marked 
"revised. " 

Sincerely, 

Hellmut H. Schmid 
Director, Geodetic Research 
and Development Laboratory 

Enclosure 

cc: l/Prof. Ivan I. Mueller, 

Ohio State University 
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July 1, 1971 
(Revises Nov. 18, 1971) 


Technical Report - Contract W-13. 321 
In accordance with contract W-13,321, dated January 19? 

1971? we are submitting the following report describing the 

contents of the magnetic tapes being supplied to the NASA 

/ 

Space Science Data Center, containing data from the BC-4 
worldwide PAGEOS network. 

, Section I 

Data Fomat of the Magnetic Tapes 

The data is furnished on 1/2" BCD (even parity) 7-track 
tapes at 556 bits/inch. Each record is 80 characters long. 

The first file on each tape contains an index to the contents 
of the tape, listing event numbers and the stations associated 
with each. The remaining files each contain data from one 
event . 

A detailed description of the data formats follows: FORTRAN 

field specifications are given for numeric dat-a fields. All un- 
specified fields are alphanumeric. 

1. Type 1 - Partially Reduced Observations 


For each event : 



• 

Event Col. 1 - 

6 


Event number 

record 

7 - 

23 


Date of observation 

24 - 

47 


Satellite name 

43 - 

53 

F6.2 

Satellite radius (m) 

54 - 


11 

No. of stations viewing event 

55 - 

61 

F7.4 

x-angle') coordinates of 

62 - 

68 

F7-.4 

> instantaneous pole 
y-angle ) relative to CIO 




(sec. of arc) 
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b. For each plate: 


Rec. 1 Col. 1 - 

6 


7 - 

30 


31 - 

33 

13 1 

35 - 

36 

12 ( 

37 - 

44' 

P8.4 / 

45 - 

48 

14 ^ 

50 - 

51 

12 > 

52 - 

59 

F8.4 / 

Rec. 2 C61. 1 - 

14 

E14.7 

15 - 

28 

E14.7 

29 - 

32 

_ . * 

33 - 

36 

14 

Rec. 3 Col. 1 - 

48 

3E16.9 

49 - 

64 

E16.9 

65 - 

80 

E16.9 

Rec. 4 Col. 1. - 

16 

E16.9 

17 - 

48 

2E16.9 

49 - 

80 

2E16.9 

Rec. 5 Col. 1 - 

48 

3E16.9 

49- 

80 

2E16.9 

Rec. 6 Col. 1 - 

16 

E16.9 

Rec. 7 Col. 1 - 

80 

4E20.13 

thru 40 



Rec. 41 Col. 1 - 

48 

3E16.9 
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Station number 
Station name 


Station latitude 


Station longitude 


0 ° 

0 ‘ 

0H 

X° 

X’ 

X" 


Station temperature (°C) 
Station pressure (mm of Hg) 
Plate number 

Number of satellite images 


Local camera orientation angles (radians) 

Angle of nonperpendicularity of plate 
coordinates (radians) 

x coordinate of principal point (m) 

y coordinate of principal point (m) 

x and y coordinate scalars (m) 


x and y coordinates (relative to princi- 
pal point) of zero distortion (m) 


3 coefficients of radial distortion 

2 coefficients of decentering dis- 
tortion 


Angle of axis of lens distortion (radians 

Upper triangular part of 16x16 covariance 
matrix associated with camera orienta- 
tion parameters (by rows) (136 terms 
total) 

Transformed camera orientation angles 
(radians) 

12 . 



For each satellite image 


(Rev. 11/17/71) 


Image Col. 1 - 8 18 

record 

. 9-36 2E14.7 


Image . number 

Adjusted x,y plate coordinates (m) 


38 - 39 12 

4l' - 42 12 


43 - 52 FLO. 6 ) observation 


54 - 55 12 


57 - 58 12 


59 - 68 FLO. 6 / observation 


2. Type 2 - Fictitious Satellite Directions 


a. For each event 


Event Col. 1 - 
record 

7 - 


Event number 


Number of stations 


Total number of fictitious 


b. For each plate 


Rec . x Col . 1 “ 


Station number 


7-30 


Station name 


31 - 34 


Plate number 


35 - 36 12 


Number of usable fictitious images 


Matrix Col. 1 - 80 4E20.13 

records 


Upper triangular part of covarianc. 
matrix associated with fictitious 


directions (radians squared). 


(Note: No. of terms ($ in the upper triangular part of the matrix 
depends on total no. (n) of fictitious points: N=n(2n+1)) 


For each fictitious point 


Record Col. 1 - 2 12 


Point number 


3 - 18 E16.9 Greenwich hour angle (radians) 


19 - 34 EI 6.9 Declination (radians) 
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Section II. 


h 

*• * 

Further Reduction of Partially Reduced Satellite Image 
Plate Coordinates 

This section of the report Qutlines the computational 

steps required to correct the partially reduced satellite 

/ 

image plate coordinates for aberrations peculiar to the 
particular comparator and camera with which they are 
associated, and to propagate the uncertainties of camera 
orientation into fully correlated image coordinate co- 
variance matrices. Formulas for transforming the corrected 
satellite image plate coordinates into directional ' coordinates 
and a discussion -of external aberrations will also be given. 

1. Symbols Used for Given Data: 

0 Nominal latitude of camera station. 

x a? y a Coordinates (angles) of earth's 

instantaneous pole relative to the CIO 
mean pole of 1900-1905. 

x,y Partially reduced satellite image plate 

/ 

coordinates . 

a L’ tu L’ K L Camera orientation angles relative to 

the local station coordinate system. 

Camera orientation angles relative to the 

coordinate system used in satellite 
triangulation adjustment (transformed 
from local system). 


14 


e 


Angle of non-perpendicularity of the 


plate coordinate system. 

x p ,y p Plate coordinates of the camera principal 

point. 

/ 

‘ C x ,Cy — X and Y plate coordinate scalars. 

A x ,Ay Plate coordinates of the point having zero 

distortion (relative to the principal 
point) . 

K^,K 2 ,K^ Coefficients of radial lens distortion 

polynomial . 

Coefficients of a function describing 

decentering lens distortion. 

An angle which relates the axis of lens 

distortion symmetry to the plate 
coordinate Y axis. 

Covariance matrix associated with camera 

i 

orientation 'parameters. 

2. Satellite Image Coordinate Corrections: 

In order to be correct, the following computations 
must be done in the specified order: 

a. Correct for non-perpendicularity of axes: 

i 

x = x + ye ; e = angle in radians 


k 4 ,k 5 

H 

[°o] 
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b. Correct for radial and decentering lens 


distortions : 

1 


* 

X = X 

* 1 

y = y 


= x ' - *p - A x 

Dy = y* - y p - * y 
a 2 = D? + D ?, 


is = K-,a 2 + K P a‘ t + K,a 6 

d 1 ^ J 


2 


-§ = k 4 + K J a 

d 


A = D x cos?>k + DySin0£ 


i) t 

x = x - 


_ji D x (2D A + d^cos0 t ) 
d X t 


(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 


" = v* - (2D.A + d 2 sin0 ) 


d 


d 


y 


( 12 ) 


, Because the distortion formulas are functions of 
undistorted distances, and because the given coordinates 
are distorted, it is necessary to iterate as. follows: 

(1) After computing (5) thru (-12) initially, set 

* 11 

X = X 

* II 

y = y 


&nd repeat (5) thru ( 12 ) • 

(2) Compare | ,A R - % n _il against a pre-set 
tolerance (0.1 micrometers is typical) and 
compare I^D n “ A D n -ll against a tolerance. 
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If both differences are less than or equal 

to the tolerances, continue with next set of 

* * 

corrections. Otherwise, set x and y equal 

t 

to the latest x" and y* 1 and repeat computations 
(5) thru (12)'. 

c. Translate coordinates to principal point 


HI.- It 

X = x - x p 

(13) 

III II 

y = y - y p 

(1*0 

Correct for scale differences 


till Ilf 

XT = X 

(15) 

• III Gy III 


<< 

II 

°l 

< II 

(16) 


At this stage, following corrections in a. thru d. 
above, the satellite image plate coordinates should be 
uncontaminated by> any known, systematic errors 
originating from the measuring or camera orientation 
process. The uncertainties of these coordinates 
attributable solely to uncertainties in the camera 
orientation parameters can be expressed by the co- 
variance matrix discussed in the next paragraph. 

3. Propagation of Camera Orientation Errors: 

Henceforth, for simplicity’s sake, let x n " and y" M 
be denoted by x and y, remembering that these symbols 
now stand for corrected satellite image plate coordinates 
and not, as originally, for given coordinates. 
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Direction Cosines 


M- 


*11 

R 12 

R 13 

R 21 

R 22 

R 23 

R 

R 

R 

31 

32 

33 J 


R 


11 
r 21 


R 


31 


R 12 

R 22 


R 


32 

R 13 

r 23 

R 33 

[kJ 


= “COs a cosK+.sinasinuusinK 
= "CosuusiriK 

” • sinacosK+cosasinuusi nn 
= -cosasinft-sinasinuucosH 

= COSUUCOSK 

= sinasinH-cosasinuucosK 
= sinacosa j 
= sinu) 

= cosacosuu 


(17) 


(18) 


Matrix computed with local orientation 
angles. 


J"Rml = Matrix computed with transformed 


orientation angles, 
b. Computational Auxiliaries 

(Note: directions cosines are from 


® = <fir 


[h] > 


c y 


(2) = R i3~ R ii 


(19) 

( 20 ) 
( 21 ) 
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L. 


© 

- (X) ^23 “^21 

( 22 ) 

© 

= © K 13- R 12 

(23) 

© 

= © R 23- R 22 

(24) 

© 

= ®R 2 1 -©R 2 2 

( 25 ) 

© 

= X-A x 

( 26 ) 

® 

= 

(27) 

© 

= ® 2 + ®' 2 

( 28 ) 

© 

= [2k 1 +4k 2 © + 6 k, © 2 ] 

(29) 

© 

= [si% +l *5 © ] . 

(30 

© 

= k x © +k 2 © 2 +k, © 3 

(3D 

© 

= k^ @ +k ? @ 2 

(32) 

© 

= k 4 +( © 2 +3® 2 )k ? 

(33) 

@ 

= (9) cos0^.- 0 sin0 t 

(34) 

© 

= 0 cos0^.+ 0 sin0 t 

(35) 

Partial derivatives 


5x 

5a 

= -C x (®@+@> 

(36) 

£z 

5a 

= -Cy( © © - © ) 

(37) 

5x 

5w 

= C x ["(l+ ® 2 )sirw L - ® ©cos©] 

(38) 

iz 

5uu 

= -C y [(l+© 2 )cos* L - ® ©sinx L ] 

(39) 

dx 

5k 

= - c x@ 

(40) 
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,10 


If = c y © 

dx 


de 


= -y 


If = 0 


dx _ i 
dx 


§* = 0 
dx 

p 

¥ = o 

By 

P 

iz = ! 

By 1 


dx 

dC 


= ® 


x 


By 

dC, 


dx 

dC 


dC 

y 

dX 

BA, 


= 0 


= 0 


2 


V (2 ® 2+ (2) 2 ) k 5 cos0 i 


£y 

BA, 


© - d; 

-4 

L 

-2 (2) @ sin0 t 

(O) - (S) 12) sin0 


COS0 


-2®{[lC4+(3 
+2 (9) kcrsin0 


2^. /r^ 2 


>*J 


COS0 

5J t 




(41) 

(42) 

(43) 

(44) 

(45) 

(46) 

(47) 

(48) 

(49) 

(50) 

(51) 


(52) 


(53) 


dX 

dA 


•IT - & (lZ cos0 t -2 (8) ( 
+'2(B) @k^cos0 t ) 


sin0 , 


- 20 - 


(54) 
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§f = -@-(D 2 @-®@ sin0 t -2 ® © cos0 t 
-4© [k^+C (B). 2 +2 ® 2 ).kjsin0 t 

(55) 

^ = ® © 


(56) 

1^=® © 

/ 

(57) 

lf 2 =® @ 2 


(58) 

ff 2 =®© 2 . 


(59) 

|g 3 =®© 3 


(60) 

li 3 =@ © 3 


(61) 

Sx = @ cos0 t +2® © 


(62) 

|f = @ sin0 t +2® © 


(63) 

M = ,fn Sx - 

Sk^ ^ Sk^ .. . ' 

; - 

(64) 

i 5 =®g 4 - ' ' 


(65) 

ff t = 


(66) 

8 t = — t] - 

- 

(67) 
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Plate coordinate covariance matrix 
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dx dx dx dx 
da du> &h de 

dy dy dy dy 
.da dm de 


dx dx dX dX dX dx dX dx dX dX dX dX 

dXp dy p dC x dC y d^ x dA y dk^ dk 2 dk 3 dk^ dk^ d0 t 

dy dy dy dy dy dy dy dy dy dy dy dy 

dx p dy p dC x dC y dA x SA y dk-j_ d ^2 dk^ d0 t 


( 68 ) 



[4 0 >[|f ] [ »o 2 ] [If ] T 

(2X2) (2x16) (16x16) (16x2) 


for a single satellite image 


(69) 


or 


2 x 

[ ax y 0 J - [ fo^ ] [ a o ] [fu^] for n satellite images 
(2nx2n) (2nxl6) (16x16) (16.x2n) 

.. "• (n < Q)' 


4. Satellite Directions and External Aberrations: 


(70) 


a 


Azimuth from south (A) and zenith distance (Z-) 

r- -i r -i 



(7D 


if U = 0 and V = 0, .A = | 

if U = 0 and V > 0, A = 3H 

* 2 

K tan (y) 

if U < 0 and V > 0, A = K + 2n 

if U < 0 and V = 0, A = K 

if U > 0, A = K + tt 


(72) 
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if W = 0, Zg = | (73) 

_ 1 _ 

if W t 0, = tan" 1 [1 s!+v£L.] 

b. Atmospheric refraction 

It is customary at this stage / to correct the zenith 
distance (Z^.) for atmospheric refraction before pro- 

X\ 

ceeding to the next step. 

c. Local hour angle (H) and declination ( 6 ) relative 

to instantaneous pole ' 

Z = Z^ corrected for atmospheric refraction (7^) 

A 

6 = sin ^ (cosZsin0-sinZcosAcos0) (75) 

Y = sinZsinA (76) 

X = cosZcos0+sinZcosAsin0 (77) 

if X = 0 and Y = 0, H = § . ( 78 ) 

if X = 0 and Y < 0. H = 3E 

2 

K = tan" 1 (I) 

if X > 0 and Y <0, H = K+ 2 n 

if X > 0 and Y > 0 , H = K 

if X < 0, H = K+tt 

d. Phase angle 

Corrections to H and 6 for phase angle, the dis- 
placement of the sun’ s reflected image from the 
satellite center, are most conveniently made at this stage. 
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e. Right ascension 

The topocentric right ascension at the time of 
observation can be computed as: local sidereal time 
of observation minus the local hour angle (H). 

5. Plate Coordinates from Corrected Directions: 

i ' 


cosZ = sin0sin6+cos0cos6cosH 
• S 


= cos0sin6-sinocos6cosH 

t 




cosZ 

_ -cos6 sinH 


m 

h 

q 


- 


cosZ 


§ 

r\ 

i 


X °Xq 


( 79 ) 

( 80 ) 

( 81 ) 

( 82 ) 

( 83 ) 

( 84 ) 


6. Greenwich Hour Angle (H 1 ) and Declination (6 1 ) Relative 
to CIO Mean Pole (1900-1905) : 


f R 0oXo] 


-cos\ 0 sino 0 

'sin\ * 
0 

C°sl 0 COSO 

-sinXo sin0 o 

-cos\ 0 

sinX Q coso 

COS0J 

0 

sino ( 

19"3861 N; \ 0 

,.0, ! »l 

= 76 4 15.1266 W 


( 85 ) 


rU**i r nr "1 • r x *i 

fc] - MW UJ 


Repeat computations (72) and (73) using U*,V* and W* 
to obtain A* and Z* , then: 

H* = A* + n (if H’ = 2ny subtract 2rr) 


6 ' = H - Z* 
2 
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2. 2 Investigations Related to the Problem of Improving Existing 
Triangulation Systems by Means of Satellite Super-Control Points 

This investigation was completed during this reporting period, the details of 
which will be published separately under Reports of the Department of Geodetic 
Science. 

The objective of this investigation was to answer the question: Whether any 
significant increment to accuracy could be transferred from a super-control net 
(continental satellite net or super-transcontinental traverse) to the basic geodetic 
net (First-order triangulation). This objective was achieved by evaluating the 
position of accuracy improvement for a triangulation station, which is near the 
middle of the investigated geodetic triangulation net. 

For the purpose of the present investigation, the triangulation of the western- 
half of the United States has been considered, as this is more accurate than that 
of the eastern-half of the United States [Simmons, 1950, p.54] . The investiga- 
tion is done on the chain from Moses Lake, Washington to Chandler, Minnesota 
(Figure 1), as these two stations are also common on both the continental 
satellite net (CSN) and the super-transcontinental traverse (STT). The data was 
supplied by the Triangulation Branch of Geodesy Division, and the Geodetic 
Research and Development Laboratory,- both of the National Oceanic and Atmo- 
spheric. Administration, Washington. 

The geodetic triangulation net is adjusted as an independent or free net, as 
it is not connected with other nets. For its unambiguous determination, besides 
the observed data which includes directions, bases (to provide the scale) and 
astronomical observations, i. e. , longitude and azimuth (to provide orientation 
of the triangulation net upon a mathematical surface, i. e. , ellipsoid), one 
fixed station is required to serve as origin [Gotthardt, 1968, p. 167]. Moses 
Lake station is kept as origin with its coordinates obtained from satellite tri- 
angulation results; these coordinates have been assumed to be the best known 
coordinates. 
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Combining the free triangulation net with super-control net of zero' order, 
i. e. , continental satellite net and/or super transcontinental traverse means 
constraining the scale and/or orientation of the triangulation net. The effect of 
this combination is comparable with "tennis racket and string effect", where 
the rigid outer racket frame (super-control) constrains the loose strings 
(triangulation net). If the strings are already constrained, there will be no 
"visible" effect of the additional constrain from the rigid outer frame. This is 
also the purpose of this investigation, i. e. , to evaluate whether the existing 
geodetic triangulation is sufficiently "constrained" or needs to be constrained by 
additional super-control net. For the present investigation a triangulation station 
Chandler, which is common to the three networks, provides constraint. 

Geodetic triangulation net can be combined with the super-control net in 
either of the two ways: 

(1) By using the actual data, i. e. , by using the actual coordinates with 
their standard errors of Chandler as obtained from CSN and STT with the 
geodetic triangulation; or, 

(2) By adding a weight constraint to Chandler with its coordinates from 
the geodetic triangulation. 

For this investigation, the first way could not be used, as the super-control 
net coordinates of Chandler station are not compatible with those obtained from 
geodetic tri angulation. As such, the second way has been preferred by using 
the actual preliminary accuracy estimates for Chandler, which are 1 part in 
385,000 and 1 part in 3 million, as obtained from CSN and STT, respectively. 
Further investigations are made by using hypothetical standard positional error 
accuracy estimates of Chandler station, which are 1:400,000; 1:500,000; 
1:600,000; 1:700,000; 1:1M; 1:1, 5M. These accuracy estimates are within the 
actual preliminary accuracy estimates of super-control nets. Thus, using 
these various accuracies of super-control net, a feeling for the accuracy limit 
of super-control net, which would be necessary to improve the investigated 
geodetic triangulation, can be obtained. 
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The Conjugate Gradient Method (Cg Method) is used for adjustment, which 
uses the original homogenized observation equations. Cg-Method has been 
programmed in such a way so as to use varying data with only change in the 
dimension statement. 

The results of the investigation are given in Tables 1 and 2, wherein the 
improvement of the particular geodetic triangulation by super-control net is 
visible only when its accuracy is at least 1 part in 500, 000. The positional 
improvement of Wyola (95), which is in tae middle of the triangulation chain, 
using various station constraints for Chandler (3) relative to free net adjustment 
are shown in Figure 2. As the preliminary accuracy of continental satellite net 
is lower than 1 part in 500, 000, this cannot be useful as a "constraint" to the 
geodetic triangulation net. On the other hand, the high accuracy of super- 
transcontinental traverse, which is one part in 3 million, makes it most 
suitable as a "constraint" to the geodetic triangulation net. 

Worth mentioning is that the longitude terms, which are Qy y anda y 2 in 
Table 1, remain practically uneffected during the entire investigation. This 
could be explained by the fact that station Wyola is very close to Laplace stations, 
which control the azimuth error accumulation, thus effecting the longitude error 
[ Bomford, 1962, pp. 90, 519]. Hence, due to closeness of Laplace stations, 
the longitude terms remain practically uneffected. 

The super-control net, i. e. , continental satellite net or super-transconti- 
nental traverse, can provide a useful constraint to the investigated geodetic 
triangulation net, and thus can improve it only when the accuracy of super- 
control net is at least 1 part in 500, 000 in this case, this corresponds to 
± 3.7m standard position error to the station Chandler. The preliminary accuracy 
of super-transcontinental traverse is already better than this limiting accuracy 
of 1 part in 500, 000. The preliminary accuracy of continental satellite net is, 
however, lower than the limiting accuracy of 1 part in 500, 000; the preliminary 
standard positional error for Chandler as obtained from continental satellite net 
corresponds to ±4. 8m, i. e. , 1 part is 385, 000. The future will show whether 


28 



Table 1 


Experiment 

Number 

Accuracy 
1 in 

■ 

H 

WYOLA (95) 

Remai'ks 

Qxx 

Qyy 

„ 2 
m x 

m y 3 

1 


2.42 

6.0 

0.5 

35.2 

2.9 

Free Net 

2 

300,000 

2.41 

6.7 

0.5 

38.9 

2.9 


3 

400,000 

2.41 

5.9 

0.5 

34.3 

2.9 


4 

500,000 

2.41 

4.1 

0.5 

23.8 

2.9 


5 

600,000 

2.41 

4.1 

0.5 

23.8 

2.9 


6 ' 

700,000 

2.41 

4.1 

0.5 

23.8 

2.9 


7 

1,000,000 

2.41 

3.7 

0.5 

21.5 

2.9 


8 

1,500,000 

2.41 

3.2 

0.5 

18.6 

2.9 


9 

3,000,000 

2.41 

2.1 

0.5 

12.2 

2.9 



Qxx> Qyy and m x 3 , m y 3 are given in 10 -4 seconds 3 . 
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Table 2 





WYOLA (95) 


Experimei 

Number 

Accuracy 




Positional Improvement 
Relative to Experiment 1 

1 in 

m x 

my 

m p 

Meters 

% 

m 

Free Net 

1.83 


1.9 



2 

300,000 



2.0 

-0.1 

- 5 

3 

400,000 

1.81 


1.8 

0.1 

5 

4 

500,000 

1.51 


1.5 

0.4 

21 

5 

600,000 

1.51 


1.5 

0.4 

21 

6 

700,000 

1.51 


1.5 

0.4 

21 

■ 

1,000,000 

1.43 


1.5 

0.4 

21 


1,500,000 

1.33 


1.4 

0.5 

26 

9 

3,000,000 

1.08 


1.1 

0.8 

42 


Standard Errors of Unknowns (m x ,m y ) and Standard Positional 
Error (m,,) are given in meters. 
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CHANDLER ACCURACY 














this limiting accuracy could be achieved by continental satellite net, especially 
because numerous spatial triangulation of CSN have produced accuracies within 
the range of 1 part in 400, 000 and 1 part in 700, 000 [Schmid, 1965, p. 22]. 
Schmid [1970, p. 23-24] indicates that continental satellite net will fall short of 
an optimum solution with respect to both its coverage and its accuracy. The 
non- achievement of the anticipated accuracy may be due to (1) the earlier 
observations made with the 300 mm lens, which have a slight disadvantage when 
compared to the presently available 450 mm lens optimized for satellite 
triangulation, (2) the lack of knowledge of the minute orientation changes of 
BC-4 system in the earlier phase, which occur occasionally between the period 
separating pre- and post-star calibrations, and (3) the absence of an optimal 
target after the demise of Echo I and II satellites, only the PAGEOS satellite is 
available, which, with an average slant distance of six million meters yields 
results degraded by a factor of three, compared with results, obtainable with a 
balloon satellite at the optimum height of 1500 km above the earth. Thus, the 
three-dimensional positions of CSN-stations will probably be determined to no 
better than ±4 meters in all components which does not seem to be good enough 
at least for this particular triangulation chain. It might be useful to have a 
"block constrain" instead of "chain constrain", that is, to use four well- 
s eparated satellite stations 003, 102, 112, and 134 (Figure 1) over a very 
large area, thus constraining the triangulation of the western-half of the United 
States instead of one triangulation chain ("chain constrain") between stations 
003 and 102. 

Super-transcontinental traverse can provide a better constraint, if more 
than two of its stations are common to the stations of geodetic triangulation net. 
Also, a ,r block constrain" as explained above, might be more useful instead of 
a "chain constrain". 
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2. 3 Scaling The SAO-69 Geometric Solution 
With C-Band Radar Data (Solution SC 11) 

SC 11 is an adjustment of the SAO optical observations combined with 
pseudo observations of chord distances derived from the adjustment of the 
C-Band network. Weighted height constraints (mean sea level + SAO 69 
undulations) were impared at most stations as a means of introducing 
valuable additional information. 

14,356 optical simultaneous observations from 28 stations were 
received from the SAO. The number and distribution of these obser- 
vations are shown in Figures 1 and 2. 

The C-Band data consisted of the results published by NASA-Wallops 
Island [1]. Upon request, they also kindly supplied us with the 
correlation matrix of their adjustment from which it was easy to compute 
the variance-covariance matrix. This adjustment was developed from 
about 2,000,000 range observations in which 466 tracks of GEOS- II 
were observed by twenty-one C-Band radars. 

Only four C-Band stations were near enough to the SAO stations 
so that the relative positions could be reliably determined through first 
order triangulation. To account for any error introduced through the 
triangulation, these ties were accounted for by imposing weighted 
relative position constraints on the adjustment. 

As the geometry of the network did not always guarantee a unique 
solution, it was also necessary to impose a few weighted direction 
constraints derived from the SAO 69 adjustment [2]. These direction 
constraints imposed were between 

Addes Ababa and Shiraz 
Addes Ababa and Natal 
Natal and Pretoria 
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Natal and Villa Dolores ' 

Natal and Arequipa 

Natal and Curacao 
Organ Pass and Rosamund 

The above data and constraints were used as input in the least 
squares adjustment program. Several adjustments were attempted, 
but the final one chosen, the SC 11, used only the two C-Band chord 
lengths (Kauai, to Merritt Island and Merritt Island to Pretoria). In 
this adjustment two iterations were carried out from which a new set 
of coordinates for the SAO Network, and the four C-Band radar sites 
were obtained. 

'T'he results of our adjustment are given in Table 1. In Table 2, 
the actual differences from the SAO 69 solution are given along with 
the standard deviations of the corresponding coordinates. For com- 
parison, also included are the standard deviations obtained by the SAO. 
Entries enclosed in a box are statistically significant at the 30 level. 
The remaining differences contain zero within their confidence intervals 
and are assumed insignificant. In those cases where the difference is 
significant, the SAO and our standard deviations are of the same order 
of magnitude. 

Tables 3 and 4 give chord-distance comparison with other adjust- 
ments on these stations. The SC 11 solution is compared with both 
SAO 69 and with G3FC 1971 [3]. 

Tables 7-8 give the residuals between the above solutions after 
systematic differences due to the various coordinate systems and scales 
were removed. Tables 5-6 contain the transformation parameters and 
the corresponding variance-covariance matrices. 

In conclusion, we can say as a result of this adjustment that the 
SAO and C-Band adjustments at least in the western hemisphere are 
compatible with one another and that the latter can successfully be used 
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to establish the scale for the former. Our standard deviations are not 
much different than those obtained by the SAO using considerably more 

observational information and the meshing together of two totally inde- 
pendent adjustments created no undue stress. 

Because of the weak Baker-Nunn connection across the Atlantic and 
the absence of a C-Band length which could be used in the eastern 
hemisphere, the results of our adjustments are weaker there than the 
SAO, but not excessively. The weakest point in the adjustment was 
Pretoria, South Africa which had optical observations only on a single 
line north to Addes Ababa and the C-Band length Merritt Island to 
Pretoria. For similar reasons, Tokyo is also poorly determined. The 
geometry of the network is not the best for geometric adjustments. 

After the coordinate transformations between the SC 11 and SAO 69 
solutions (Tables 5-6), it seems that after the removal of the expected 
differences in the coordinate systems a small scale difference (1 ± .17 ppm) 
is evident in the global network, which, considering that the SC 11 scale 
is based on radar measurements and the SAO 69 on the adopted value of 
l^M, is a very satisfactory agreement. On the other hand, in the 
European part of the network a large scale difference (-11 ± 2 ppm) is 
evident. None of the residuals in Tables 7-8 are statistically significant 
on the 3 a level, those on the 2 a level are boxed. 

Further details on this and other solutions may be found in a separate 
report to be published soon. 
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Figure 1 

SAO 69 Simultaneous Global Observations 




Figure 1 (continued) 
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Table 1 (continued) 
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Table 2 

Coordinate Differences SC 11 - SAO 69 



All units are meters. 
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Table 4 

Differences Between Distances GSFC 71 - SC 11 
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All distances in meters 
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Table 5 

Global Transformations SC 11 - SAO 69 
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All rotations are positive for counterclockwise rotation as 
viewed looking toward the origin from the positive end of 
the rotation axes as follows: 
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Table 6 

European Transformations SC 11 - SAO 69 
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Table 7 

Global Residuals SC 11 - SAO 69 



Transf. 3 Param. ! 

Transf. 7 Param. 

Station 

AX 

AY 

AZ 

AX 

AY 

AZ 

9001 

- 8.1 

9.1 

- 1.6 

-10.2 

5.0 

- 1.8 

9004 

- 6.6 

3.7 

4.8 

- 4.6 

8.2 

8.9 

9006 

7.4 

- 1.7 

4.6 

2.2 

6.7 

4.4 

9007 

-11.2 

3.2 

-17.4 

- 5.8 

- 0.3 

-20.5 

9008 

- 3.4 

0.5 

4.0 

- 5.6 

9.1 

5.5 

9009 

- 6.6 

- 2.6 

4.9 

- 2.9 

- 5.5 

5.2 

9010 

- 2.9 

- 5.7 

4.2 

- 1.7 

- 8.9 

5.2 

9011 

-14.6 

- 3.5 

-20.9 

- 8.2 

- 6.1 

-25.6 

9012 

30.7 

- 0.9 

- 9.8 

23.7 

- 4.5 

-14.2 

9028 

23.3 

8.9 

- 6.9 

24.5 | 

17.6 

- 6.5 

9029 

-19.4 

1 

-12.0 

- 5.7 

-12.6 

-11.3 

- 5.8 

9031 

-15.2 

5.1 

8.4 . 

1 

- 9.1 

2.8 

2.0 

9091 

-16.0 

-22.5 

- 5.1 

-16.1 

-15.6 

- 1.6 

9113 

20.9 

13.0 

19.6 

17.4 

9.0 

19.0 

9114 

- 2.2 

17.2 

11.5 

- 6.2 

15.3 

13.1 

9115 

- 6.1 

5.6 

4.2 

- 7.9 

10.5 

8.5 

9117 

41.1 1 

-31.8 

1.8 

33.1 

-34.4 

- 3.6 


All distances in meters. 
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Table 8 

European Residuals SC 11 - SAO 69 



Transf. 3 Pa ram. I 

Transf. 7 Param. 

Station 

AX 

AY 

AZ 

AX 



AY 

AZ 

8015 

5.3 

- 3.1 

6.3 

5.1 

- 5.1 

- 2.4 

8019 

4.6 

- 2.7 


2.9 

- 3.5 

1.4 

9065 

-17.6 

- 6.6 


-10.7 

- 9.8 


9066 

+ 1.3 

- 8.5 


2.3 

- 9.1 

7.1 

9074 


- 4.2 


0.5 

1.7 


9077 


- 8.4 

- 6.8 

-11.3 

- 0.0 

- 9.4 

9080 

- 8.3 

40.5 | 

-30.9 

2.7 

32.6 

-12.4 


All distances in meters 


















r 


2.4 Geodetic Satellite Observations in 
North America Solution NA-8 

The results of the adjustments of the GEOS-I tracking network through 
the NA-6 adjustment were reported in [1]. After the completion of these 
adjustments, it was apparent that the weakest of the adjusted station coor- 
dinates were the heights. The approximate height coordinates used in the 
adjustments were taken from the station descriptions on the geodetic data 
sheets of [2]. The only height constraint imposed was at Columbia, 
Missouri; all others were allowed to adjust freely. At the time, it was 
not possible to constrain any station height with any degree of accuracy. 

After the completion of all previous adjustments, a new geoid became 
available from SAO [3] . This geoid gave the heights above the SAO 
ellipsoid. This ellipsoid is earth-centered, and based on comparison of 
station coordinates in the continental United States, the following shifts 
were determined for the North American Datum: 

Ax = - 38 m 

Ay = 164 m 1 

Az = 175 m 

The sign convention of these shifts is SAO-NAD. 

With the geoid map, it was possible to determine the geoid undulations 
at each of the observing stations in the optical network. Since the ortho- 
metric heights were well determined at the stations, it was simply a matter 
of adding the geoid height and the orthometric height to arrive at heights 
with respect to the SAO ellipsoid. By performing a datum transformation 
the heights were computed with respect to the Clark 1866 ellipsoid, i.e. , 
with respect to the NAD. 

The NA-6 solution was readjusted with the new computed heights as 
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constraints placed on all 30 optical stations, using a weight corresponding 
to a standard deviation of 5 meters. This was referred to as the NA-8 
solution; the results are listed in Table 1. 

The NA-8 solution shows the adjusted coordinates to be realistic, 
the standard deviations of the adjusted coordinates being smaller than those 
of any other OSU adjustment. The NA-8 coordinates cannot be compared 
directly with the NAD coordinates because of the height change at the 
origin. In order to make a comparison with the NAD coordinates, the 
following coordinate differences must be added to the NA-8 coordinates: 

Ax = - 1 . 6 meters 
Ay = +20.4 meters 
Az = -20. 5 meters 

These are the shifts of Columbia, Missouri from its NAD coordinates. 

The NAD geodetic coordinates of the stations are listed in Table 2. 

Table 3 shows the datum transformation parameters NA8-NAD. The 
interpretation of the parameters are identical to those in [1]. Compari- 
sons were also made with other solutions, namely, those of the Smith- 
sonian Astrophysical Observatory [3], the Goddard Space Flight Center [4], 
and our previous NA-6 [1], 

The comparisons show that after removing the systematic differences 
due to the different coordinate systans and scale (Table 4), the residuals 
are smaller than those expected from the standard deviations on the 3 a level 
(Table 5). Those residuals which are greater than 2a (but still smaller than 
3a) are framed in Table 5. 2a in this case corresponds approximately to 
a confidence of 95%. 

More details and other computations will be published in a separate report. 
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Coordinates of the North American GEOS-I Tracking Stations from the NA-S Geometric Adjustment 
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All coordinates and standard deviations in meters. 
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General Information: No. of degrees of freedom 5254 

No. of ground stations 34 Quadratic sum of the residuals (V'PV) 5082 

No. of spatial chord equations 1 Standard deviation of unit weight 1 
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Note: The above coordinates were arrived at by applying the shifts AX = -1.6m, AY - 29.4m and 

AZ = -20.5m to the NA-8 coordinates and then converting these values to ellipsoidal coor- 
dinates on the ellipsoid a = 6378206.4, b = 6356583.8. 
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V 


Table 3 

Datum Transformation Parameters: NA-8-NAD 



Eastern Half 
(14 stations)* 

Western Half 
(5 stations)** 

Combination 
(19 stations) 

Veis 

M") 
(") 
e x C) 

-1.2 ± 0.4 
-0.2 ± 0.4 
1. 6 ± 0.5 

0.8 ± 0.8 
1.0 ± 0.9 
-0.6 ± 1.0 

-1.0 ± 0.3 
-0.1 ± 0.4 

0.6 ± 0.2 
« 

Molo- 

densky 

e * (”) 

0y o 

e*n 

-2.0 ± 0.5 
0.0 ± 0.4 
-0.2 ± 0.4 

0.9 ± 1.0 
-0.4 ± 1.0 
1.0 ± 0.9 

-1.1 ± 0.3 
0.4 ± 0.3 
-0.1 ± 0.4 


e (x 10" 5 ) 

1.3 ±1.9 

-0.4 ± 4.2 

0.1 ± 1.4 

Shifts 

AX(m) 
A Y (m) 
A Z(m) 

1.0 ± 2.9 
-31. 1 ± 3.2 
20.8 ±3.1 

8.2 ± 4.6 
-28.8 ± 4.5 
12.6 ±4.2 

3.1 ± 2.0 
-25.3 ± 1.7 
17.0 ± 1.7 


♦Eastern Stations: 
♦♦Western Stations: 


1021, 1022, 1034, 1042, 3334, 3401, 3402, 
3648, 3657, 3861, 3037, 7043, 7072, 7075. 

1030, 3400, 3902, 7036, 7045. 


All rotations are about Meades Ranch and are positive for 
counter-clockwise rotation as viewed looking toward the origin 
from the positive end of the rotation axes as follows: 


Veis: « 


Molodensky : 


x horizon plane south 
y horizon plane east 
z zenith 

x parallel to 0 meridian 
y parallel to 90° E meridian 
z parallel to rotation axis of 
the earth (ICO) 
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Transformations Between Various Solutions 
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In these transformations, the shifts mentioned on page 51 have already 
been taken into account, i.e., the correct heading should be for example 
(NA 8 NAD) - GSFC, etc. 




Station- coordinate Residuals After Transformations 
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3. PERSONNEL 


Ivan I. Mueller, Project Supervisor, part time 

Festus Charles, Graduate Research Associate, part time (September only) 

Muneendra Kumar, Graduate Research Associate, part time 

James P. Reilly, Graduate Research Associate, part time 
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Marvin C. Whiting, Graduate Research Associate, part time 

Christine DeGeorge, Research Aide, part time 
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4. TRAVEL 

Ivan i. Mueller 

Key Biscayne Florida, October 5-8, 1971 
To attend NOAA/NASA/NAVY Conference . 

Ivan I. Mueller 

Greenbelt, Maryland, November 17-19, 1971 

Discussions with personnel at NASA, Goddard Space Flight Center 
regarding project reporting. 

Ivan I. Mueller 

San Francisco, California, December 5-7, 1971 

Present two papers at the Annual Fall Meeting of the American 

Geophysical Union. 
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by Hans D. Preuss 
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71 Proposed Optical Network for the National Geodetic Satellite Program 
by Ivan I. Mueller 

May, 1966 

82 Preprocessing Optical Satellite Observations 
by Frank D. Hotter 
April, 1967 

86 Least Squares Adjustment of Satellite Observations for Simultaneous 
Directions or Ranges, Part 1 of 3: Formulation of Equations 

by Edward J. Krakiwsky and Allen J. Pope 
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87 Least Squares Adjustment of Satellite Observations for Simultaneous 
Directions or Ranges, Part 2 of 3: Computer Programs 

by Edward J. Krakiwsky, George Blaha, Jack M. Ferrier 
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88 Least Squares Adjustment of Satellite Observations for Simultaneous 
Directions or Ranges, Part 3 of 3: Subroutines 

by Edward J. Krakiwsky, Jack Ferrier, James P. Reilly 
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93 Data Analysis in Connection with the National Geodetic Satellite Program 
by Ivan I. Mueller 
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No. NGR 36-008-093: 

100 Preprocessing Electronic Satellite Observations 
by Joseph Gross 
March, 1968 

106 Comparison of Astrometric and Photogrammetric Plate Reduction Techniques 
for a Wild BC-4 Camera 
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110 Investigations into the Utilization of Passive Satellite Observational Data 
by Janies P. Veach 
June, 1968 

114 Sequential Least Squares Adjustment of Satellite Triangulation and 
Trilateration in Combination with Terrestrial Data 
by Edward J. Krakiwsky 
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118 The Use of Short Arc Orbital Constraints in the Adjustment of Geodetic 
Satellite Data 
by Charles R. Schwarz 
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125 The North American Datum in View of GEOS I Observations 
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139 Analysis of Latitude Observations for Crustal Movements 
byM.G. Arur 

June, 1970 

140 SECOR Observations in the Pacific 

by Ivan I, Mueller, James P. Reilly, Charles R. Schwarz, Georges Blaha 
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147 Gravity Field Refinement by Satellite to Satellite Doppler Tracking 
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148 Inner Adjustment Constraints with Emphasis on Range Observations 
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150 Investigations of Critical Configurations for Fundamental Range Networks 
by Georges Blaha 
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The following papers were presented at various professional meetings: 
"Report on OSU participation in the NGSP" 

47th Annual meeting of the AGU, Washington, D. C. , April 1966 
"Preprocessing Optical Satellite Observational Data" 

3rd Meeting of the Western European Satellite Subcommission, IAG, Venice, 
Italy, May 1967. 

"Global Satellite Triangulation and Trilateration" 

XIVth General Assembly of the IUGG, Lucerne, Switzerland, September 1967, 
(Bulletin Geodes ique , March 1968). 

"Investigations in Connection with the Geometric Analysis of Geodetic Satellite 
Data" 

GEOS Program Review Meeting, Washington, D. C. , Dec. 1967. 

"Comparison of Photogrammetric and Astrometric Data Reduction Results for 
the Wild BC-4 Camera" 

Conference on Photographic Astrometric Technique, Tampa, Fla. , March 1968. 

"Geodetic Utilization of Satellite Photography" 

7th National Fall Meeting, AGU, San Francisco, Cal., Dec. 1968. 

"Analyzing Passive- Satellite Photography for Geodetic Applications" 

4th Meeting of the Western European Satellite Subcommission, IAG, Paris, 

Feb. 1969. 

"Sequential Least Squares Adjustment of Satellite Trilateration" 

50th Annual Meeting of the AGU, Washington, D. C. , April 1969. 

"The North American Datum in View of GEOS-I Observations" 

8th National Fall Meeting of the AGU, San Francisco, Cal. , Dec. 1969 and 
GEOS-2 Review Meeting, Greenbelt, Md. , June 1970 (Bulletin Geodesique, 
June 1970). 

"Experiments with SECOR Observations on GEOS-I" 

GEOS-2 Review Meeting, Greenbelt, Md. , June 1970. 

"Experiments with Wild BC-4 Photographic Plates" 

GEOS-2 Review Meeting, Greenbelt, Md. , June 1970. 

"Experiments with the Use of Orbital Constraints in the Case of Satellite Trails 
on Wild BC-4 Photographic Plates" 

GEOS-2 Review Meeting, Greenbelt, Md. , June 1970. 
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"GEOS-I SECOR Observations in the Pacific (Solution SP-7)" 

National Fall Meeting of the American Geophysical Union, San Francisco, 
California, December 7-10, 1970. 

’’Investigations of Critical Configurations for Fundamental Range Networks" 
Symposium on the Use of Artificial Satellites for Geodesy, Washington, D. C., 

April 15-17, 1971. 

"Gravity Field Refinement by Satellite to Satellite Doppler Tracking" 

Symposium on the Use of Artificial Satellites for Geodesy, Washington, D. C. , 

April 15-17, 1971. 

"GEOS-I SECOR Observations in the Pacific (Solution SP-7)" 

Symposium on the Use of Artificial Satellites for Geodesy, Washington, D. C. , 

April 15-17, 1971. 

"Separating the Secular Motion of the Pole from Continental Drift - Where and 
What to Observe?" 

IAU Symposium No. 48, "Rotation of the Earth, ” Morioka, Japan, May 9-15, 1971. 

"Geodetic Satellite Observations in North America (Solution NA-8)" 

Annual Fall Meeting of the American Geophysical Union, San Francisco, 
California, December 6-9, 1971. 

"Scaling the SAO-69 Geometric Solution with C-Band Radar Data (Solution SC 11)" 
Annual Fall Meeting of the American Geophysical Union, San Francisco, 
California, December 6-9, 1971 
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